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Abstract 

We propose the new mechanism of neutrino flavour relaxation to 
explain the experimentally observed changes of initial neutrino flavour 
fluxes. The test of neutrino relaxation hypothesis is presented, using 
the data of modern reactor, solar and accelerator experiments. The fi- 
nal choice between the standard neutrino oscillations and the proposed 
neutrino flavour relaxation model can be done in future experiments. 
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Now the phenomena of changes in initial neutrino flavour flux is observed 
in different neutrino experiments. The SNO experiment pQ evidently detects 
only 1/3 part of initial electron neutrino flux from the Sun and 2/3 part 
of muon and(or) tau neutrino fluxes. The Kamland reactor experiment [2] 
detected only 61% of expected electron antineutrino events from different 
reactors at mean distance of 180 km. Convincing evidence of initial neu- 
trino flavour flux changes is observed also in Super-K j3] and MACRO 0] 
atmospheric neutrino data and K2K j3] accelerator experiment with muon 
neutrinos. 

The standard common way of interpreting these results lies in neutrino 
oscillation hypothesis, first proposed by Bruno Pontecorvo jU] and developed 
in further works [7j. 

Here we discuss the alternative mechanism of neutrino flavour relaxation, 
which can also describe the observed changes of initial neutrino flavour fluxes 
with distance. 

The proposed model is similar to the mechanism [S] of spin relaxation 
in random fluctuating magnetic field B with zero average < B(t) >= 
and mean square fluctuating field value < B 2 (t) >^ 0. The spin relaxation 
process is described by W.Pauli master equation 0. 

Let us assume the existence of some small random fluctuating vacuum 
field V, causing the transitions between different lepton flavours. Such field 
V can have mean zero value, with < V 2 (t) > different from and (< V 2 > 
) x / 2 > m Ue , m v , m Vr . Interactions of neutrino flavours with the vacuum field 
should lead to flavour relaxation process. 

The time evolution of the neutrino states are governed by the Shrodinger 
equation: 



where | u(t) > is the neutrino vector of state and H(t) is the time- 
dependent Hamiltonian of the system, which form depends on in what basis 
it is given. 

In flavour basis the total Hamiltonian in the random fluctuating vacuum 
field for rest reference frame can be written as 



i-\v(t) >=H(t)\u(t)>, 



(1) 



H f (t) = H m + V f (t) 



(2) 
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where 
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(3) 

Here H m is the free Hamiltonian in mass basis. Note, that in this model 
the neutrino flavour states are assumed to be the same as their massive states 
(or in other words the mixing mass matrix is diagonal). 

V aa i (a = e, fi, t) is the vacuum field potential with mean value < V aa >(t) > 
and < {V aa ,{t)f >^ 0. 

The neutrino flavour evolution in time can be obtained from Eq. (1-3), 
using the density matrix approach ^U] and is given by Pauli master equation: 
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where n a (t) (a = u e , u^, u T ) are the probability of observing neutrino with 
electron, muon or tau flavour, W aa < corresponds to the neutrino transition 
rates from flavour a' to flavour a and 



E n a (t) = l 



a=v e ,vu,,v r 



The general solution of Eq. (4) is given by the sum of two exponents and 
constant: 



t/Tl + c 2 e-W , (5) 






with 

3 3 3 

2^ = 1, £> = 0, ^ Cj = 0. (6) 

i=l i=l i=l 

Note, that in relaxation model the neutrino flavour lepton numbers are 
violated, while the total sum of the numbers are constant (also as in the case 
of oscillation model). 
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To illustrate the possibilities of the relaxation model the particular " sim- 
ple" solution can be found under the following assumptions: 

a) = W eT = W Te , W^ T = W T ^ 

b) > W efl , W eT 

Now equation (jSJ) looks as: 



/ n v M \ ( 1/3 \ / \ / -1 \ 

\n v {t) \= 1/3 + 6 1 e- 2W ^ + c\ 1/2 e -m^w„)t 

\nJ T (t)J V 1/3 / V- 1 / V 1/2 / 

(7) 

For small mass differences between different neutrino flavours {m Ve ~ 
m v vh Vt and Am u <C m v ) with 7 = E u /m v ^> 1 the probability P of 
observing the neutrino flavour with energy E u at distance r from the source 
is given by: 



PuAr) \ ( 1/3 \ / \ / -1 \ 

P^(r) = 1/3 +6 1 e -'/<Ao^) + c 1/2 e-rMiK). 

^(r)/ \ 1/3 / V- 1 / V 1 / 2 / 

(8) 

where we use the new notation A = (2W flT m u ) 1 
and A x = 2/3[(W e ^ + W eT )m v \-\ 

Equation (|SJ) can be solved for different initial experimental conditions of 
neutrino flavour fluxes. 

For initial pure electron neutrino or antineutrino flux (like in the case of 
solar or reactor experiments) P Ue (0) = 1, P v A0) = Pv T (®) — the solution is 
given by: 

P Ue (r) = l/3[l + 2c- r /^)], 
P„» = 1/3[1- e-/(^)], 

P„ T (r) = 1/3[1- e^M^)]. (9) 

From Kamland experiment data [2] it is possible to estimate the value 
of Ai parameter. Taking the effective energy of reactor antineutrino flux 
to be equal to 4.8 MeV (here we take into account the threshold of Kam- 
land detector Ey e > 3.4 MeV), mean reactor distance 180 km and R = 

Pu e (measured) / Pv e (expected) = 0.61 We obtain 

A l = 43 km I MeV . 
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Fig.l illustrates the distance dependence of reactor antineutrino fluxes ( 
for = 4.8 MeV). 
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Figure 1: Probability to observe different flavours from reactor for Ep = 
4.8 MeV. 



Note, that the final flavour survival probability 1/3 is consistent with the 
SNO [Tj experimental data with the ratio of measured neutrino flavour fluxes 

CC 

■— = 0.306 ± 0.026(stat) ± 0.024(syst) , (10) 

iv (_> 

and with the absence of distortion in the measured neutrino spectrum at 
low energies, while such distortion is predicted by LMA solution [TT] . 

It agrees also with the Homestake results for the ratio R of observed 
and predicted by SSM [T3] neutrino rates: R = 0.34 ± 0.03. 

The deviation from "simple" relaxation model exists for Gallium exper- 
iments JH] with R = 0.553 ± 0.034. But it is worth to mention, that 
the R value depends on the accuracy of Standard Solar Model predictions 
for low energy fluxes. 

For the case of pure initial muon neutrino or antineutrino flux (like in 
accelerator experiments) P u ^) = 1, Pv e {fy = ^W(O) = the solution of Eq. 
(9) is: 
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P Ue (r) = 1/3[1- e"^ Al ^)], 

p^( r ) = 1/3 + 1/2 e- r/(AA) + l/6 e- r/(AA) , 

P VT ( r ) = 1/3-1/2 e" r/(AA) + l/6 e" r/(AA) . (11) 

Preliminary estimation of Ao parameter can be done using K2K accel- 
erator experiment data 5J. Taking the mean energy of muon neutrino flux 
E v = 1.3 GeV, distance r = 250 km, R = ( meaS ured)/^ (expected) = 0.70 we 
obtain: 

A = 0.21 km/MeV . 

Fig. 2 illustrates the distance dependence of neutrino fluxes from muon 
neutrino beam [E v = 1.3 GeV). 
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Figure 2: Probability to observe different flavours from accelerator muon 
neutrino beam for E u — 1.3 GeV. 

It is necessary to note that the estimations of neutrino relaxation param- 
eters Ao and Ai are very preliminary and were made mainly for illustration 
of relaxation model. More precise calculations can be done in future, using 
exact data of different experiments (taking also into account atmospheric 
neutrino data). 
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Finally we would like to emphasize the main difference between the stan- 
dard oscillation theory and the proposed here relaxation model: the depen- 
dence of neutrino flavour fluxes from distance is described by the sum of 
constant and 2 exponents, instead of oscillation case. Fig. 3 illustrates this 
difference for reactor antineutrino experiments. 
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Figure 3: The ratio of measured to expected v e flux from Kamland j2] and 
Chooz fB] reactor experiments. Dotted curve — predictions of oscillation 
model with sin 2 29 = 0.91 and Sm 2 = 6.9 • 10 _5 (eV^) 2 best fit parameters from 
Solid curve — predictions of "simple" relaxation model. 

Future reactor and accelerator neutrino experiments can provide the nec- 
essary data to choose between the neutrino oscillations and the proposed 
flavour relaxation model. Of course, the possibility to have a mixture of 
neutrino oscillation and relaxation models also exists. 

Submitted to JETP Letters. 
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